Appearance of classical Mixmaster Universe 
from the No-Boundary Quantum State 
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We investigate the appearance of the classical anisotropic universe from the no-boundary quan- 
tum state according to the prescription proposed by Hartle, Hawking and Hertog. Our model is 
homogeneous, anisotropic, closed universes with a minimally coupled scalar field and cosmological 
constant. We found that there are an ensemble of classical Lorentzian histories with anisotropics 
and experience inflationary expansion at late time, and the probability of histories with anisotropics 
are lower than isotropic histories. Thus the no-boundary condition may be able to explain the emer- 
gence of our universe. If the classical late time histories are extended back, some become singular 
by the existence of initial anisotropies with large accelerations. However we do not find any chaotic 
behavior of anisotropies near the initial singularity. 



I. INTRODUCTION 



Quantum cosmology is one of the most ambitious sub- 
ject in physics. Its aim is to treat the universe itself by 
quantizing spacetime and to explain the creation of uni- 
verse. We have not been able to treat the subject fully 
satisfactory because we do not have a complete theory 
of quantum gravity. However we do have very interest- 
ing approaches to the subject using best knowledge of 
our understanding of general relativity and quantum field 
theory. One is the quantum cosmology based on canon- 
ical quantization of general relativity. By using Hamil- 
tonian formulation of general relativity and regarding 3- 
metric and its conjugate as fundamental operators, one 
can obtain Schrodinger like equation on the space of 3- 
metric(superspace), known as Wheeler- DeWitt equation 
H'^{q^) = 0, where the are a set of coordinates for 
minisuperspace and ^{q"^) is interpreted as wave function 
of the universe. By restricting the spacetime to homo- 
geneous and isotropic(called minisuperspace), Wheeler- 
DeWitte equation greatly simplifies and an interesting 
idea for the origin of the universe have been proposed fH. 
However, nobody knows how to calculate the probabil- 
ity of specific universe in this approach, and there are 
some problems e.g. the problem of time. There has been 
some progress in the past 10 years or so in this approach 
based on loop quantum gravity(e.g. Ashtekar Bo- 
jowald i). 

The another approach is based on the path integral 
quantization and is called no-boundary proposal for the 
wave function of the universe 4] . In this approach the 
wave function of the universe is defined by summing all 
Euclidean compact geometry weighted by exponential of 
action exp{—I) without boundary. As Hawking said ; 
the boundary condition of the universe are that has no 
boundary (Hartle & Hawking [3], Hawking 0]). Accord- 



ing to this proposal the origin of the quasiclassical realm 
from the no-boundary proposal for the quantum state of 
the universe has been discussed by Hartle, Hawking and 
Hertog in a class of minisuperspace models with ho- 
mogeneous, isotropic closed spacetime geometries. They 
argued that a wave function with the semiclassical form 



*(g'4) = A{q^)e'^'^'i^^''' 
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in some region of minisuperspace predicts an ensemble 
of classical trajectories provided that S{q^)/h is rapidly 
varying and A{q^) is slowly varying. It was shown 
that the complex solutions (metrics) of Euclidean Ein- 
stein equations that extremize the Euclidean action are 
distinct from the real Lorentzian classical metrics for 
which they provide the probabilities, and the probabili- 
ties for classical histories derived from the no-boundary 
wave function do define a measure on classical phase 
space. 

It will be very interesting to see if it is possible to ap- 
ply their method to more general class of minisuperspace. 
This is the aim of the paper. Namely we investigate the 
emergence of classical universes with anisotropy accord- 
ing to the prescription proposed by HHH. We consider 
the Bianchi type IX spacetime which is homogeneous but 
not isotropic, with a minimally coupled scalar filed mov- 
ing in a quadratic potential and cosmological constant. 
Quantum cosmology of Bianchi type IX model has been 
studied using Wheeler-DeWitte equation. For example, 
Battisti et al. 0] investigated the situation where the 
isotropic scale factor is considered as intrinsically differ- 
ent variable from anisotropies. Namely, they regarded 
the scale factor as time-parameter, and showed that the 
wave function of the universe is spread over all values of 
anisotropy near the cosmological singularity. However, 
it is also shown that the wave function asymptotically 
peaked around the isotropic configuration when the uni- 
verse grows large. We will also find that the universe 
approaches isotropic configuration, but possible range 
of anisotropy near the bounce or singularity is very re- 
stricted unlike their result, but is consistent with the pre- 
diction by loop quantum cosmology [l^. 
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The organization of the paper is as follows. We briefly 
review mixmaster universe in Sec. II and explain the con- 
dition of classicality according to HHH in Sec. III. Then 
we calculate the wave function of our model and applies 
the prescription of HHH in Sec. IV. The numerical re- 
sults are shown in Sec. V. Finally summary and some 
discussion will be given in Sec. VI. 



II. MIXMASTER UNIVERSE 

In this section we briefly review Bianchi type IX model. 
The geometry of Bianchi type IX model may be written 
as 



The conjugate momenta of a and P± are as follows. 



(2) 
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It is possible to think the Hamiltonian ([6]) as that of 
a particle with triangular shaped potential 
whose walls are very steep. Therefore the evolution of the 
closed universe with anisotropy may be expressed by the 
motion of this point mass (e.g. Misner,Thorne & Wheeler 
_8]). Since the motion of the point shows chaotic behav- 
ior in general by bouncing between potential walls, it is 
argued that the classical universe is predicted to behave 
chaotically taking all possible values of anisotropics. This 
model of the universe is called the mixmaster universe . 



where TV is the lapse function and Pij represent spacetime 
anisotropy. The 1-forms tr' are defined by 



— cos ijjd9 — cos ip sin 

— sin ^jjd9 — cos ■0 sin 9dcj) 
= dip + cos Odip 



(3) 



N, a, [3ij are functions of time t. The anisotropy variables 
(3± are parametrized as follows; 



= /3+-I-V3/3- 
/333 = -2/3+ 



(4) 



Thus the dynamical variables of our model are a.^(3±. 
This model is reduced to closed FRW model when /3± — 
0. The (Lorentzian) action of this model is calculated to 
be 



5 = - / dtNe'°' 



Ar2 V 



(5) 

where we do not include cosmological constant at the 
moment, the three-curvature "^-Rix is 



III. NO BOUNDARY WAVE FUNCTION OF 
THE UNIVERSE AND SEMI-CLASSICAL 
APPROXIMATION 

In this section we review the semi-classical approxima- 
tion of the no-boundary wave function of the universe 
according to HHH. 

The no-boundary wave function is defined by sum-over- 
histories 



'I'[q^] = J Sq'^cxp{-I[q'' 



ix)]/n) 



(9) 



where q'^{X) are histories of dynamical variables such as 
a,P±,etc, and I[q^{X)] is their Euclidean action. The 
sum is carried out over cosmological geometries that are 
regular on a manifold with only one boundary at which 
q"^(A) take some prescribed real values. The integration 
is carried out along a suitable complex contour which 
ensures the convergence of the integral and the reality of 
the result. 

We use the steepest descent approximation to evaluate 
the path integral. We consider only paths that extremize 
the action. 



=e-2"[2e-2/3+ cosh2\/3/3- 



-e-^^+ + e*f^+{l - cosh4\/3/3_) 



then we can write 
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V(A) 
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(10) 



e-^"l/(/3+,/3_) 



(11) 



Then the Hamiltonian for the Bianchi type IX universe 
is 



GN 



-3a 



r 2,2,2 

[-Pa+P++P- 



2 ' ^^^'e4«F(/3+,/3_)] 
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(6) 



where the summation run over the extrema. The equa- 
tions given by (fTO|) are recognized as Euclidean Einstein 
equations. Note that these extremized paths are com- 
plex contours, not real with a few exceptions (Halliwell 
and Hartle [9]). Of course the physical quantities are real. 
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thus we will set the boundary value as real. Therefore 
the extremized paths are fuzzy instanton, which started 
as complex Euclidean geometry and growing to asymp- 
totically real Lorentzian at large scale. 
We consider only the lowest order in h, thus 

Axt = IcAq"^) + 0{h) 

EE I[q^^,iX)]+Oih) (12) 

From now we drop the subscript 'ext'. Since the Eu- 
clidean action is complex, we can divide the action into 
real part and imaginary part, I = In — iS and write the 
wave function as follows; 

A{q^) = e-^^(«^)/'^ (13) 

The wave function is said to be the approximate semi- 
classical form when S{q^)/h varies rapidly and A{q^) 
varies slowly. Then the wave function in this region of 
minisuperspace predicts an ensemble of classical trajec- 
tories. Therefore 



STEP2 

Solve these equations in complex plane under the condi- 
tion that the boundary values are real. 

STEPS 

Check whether the solutions satisfy the classicality con- 
dition or not. 

STEP4 

Solve the classical Lorentzian paths and calculate the rel- 
ative probabilities of these histories. 

IV. NO-BOUNDARY WAVE FUNCTION FOR 
MIXMASTER UNIVERSES 

In this section we calculate the steepest descent ap- 
proximation to the no-boundary wave function for the 
mixmaster universes according to the above prescrip- 
tions. We include the cosmological constant and a mini- 
mally coupled single scalar field with quadratic potential. 

By replacing t iS ^ —I in the Euclidean 

action becomes 



\^aIr\ « \VaS\ 



(14) 



is called the classicality condition. Using conserved cur- 
rent 



(15) 



we assume that the relative probability density of clas- 
sical history is the component of the conserved current 
along the specific surface in minisuperspace. 



p{q^) EE J-n^\A{q^)\^V^S{q^) 

= e-2^«(«-')V„5 (16) 

n is the normal of the spacelike surface in minisuper- 
space. If we get the complex paths which extrcmizc the 
Euclidean action and satisfy the classicality condition, we 
are able to calculate the classical Lorentzian history by 
using these paths as initial condition. 



Lorentz {Ql 



Im[/(^^ 



pi = -lm{p^) 



(17) 



Summarize this section 
STEPl 

Calculate equation (10) to get the Euclidean Einstein 
equations. 



/KA),$(A),/3±(A)] = /ix + /$ 
where Iix is the action of spacetime 



(18) 
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with A' ~ dA/dX. /$ is the action for the minimally 
coupled scalar field with a quadratic potential 



/$ = d X^J~^Lr, 



(20) 



We have assumed that scalar field is homogeneous, and 
m is the mass of scalar field. For later convenience, we 
rescale the variables as follows; 



1- 



47rG' 
A 



A* 



(21) 
(22) 

(23) 

(24) 

(25) 



Then the explicit form of the Euclidean action takes the 
following form 
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I[a{X),^\),f3±{X)] 



9tt 
4AG 



dXNa^ 



1 

iV2 



+ (3'^ + f3'3 



+i--^{i-vm) + 



7V2 



2 J,2 



(26) 



where 



+ ie-«^+-^e4''+(l-cosh4V3/3„) (27) 

Note that the variables are all complex valued. Therefore 
the parameter A can be chosen to be real, and we choose 
A = at the starting point and A = 1 is our boundary. 

Taking the variation to the above Euclidean action and 
using the new complex parameter dr = NdX, we get five 
equations 



we can Taylor expand the variables near the origin t = 
and get the following expressions 



a(r)-r-i(l + MVo')r3 + 0(T^) 
6 

0(t) = 00 + -mVot 
1 



(33) 



/3+(r) = P2T 



m'0o(2 + + 2fi^o^y + 0{t') (34) 



^(1 + m'4')^2 - 10P2' + 10M2' 



+ 0(r«) 



/3_(r) = A#2t2 
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(1 + fi^o^)M2 + 5P2M2 



+ 0(r^ 



(35) 



(36) 



This shows that we can choose three complex numbers 
(j)o,P2,M2 as our free free parameters. We write them 
as follows: 



+ - (1 - T/(/3±)) +a2,/,2 +aV^0^ = (28) 

«^ 3 , • 2 • 2, 

a+ - + -«(/?+ +/?-) 

-^(l-n/3±)) + ^a(<^'+^V) = (29) 



+ 3-(?!)-Ai> = O 

P± + 3-/^± - 

a 2a"^ op± 



(30) 
(31) 



where A — dA/dr. Four in the above five equations are 
independent. Finally, the Euclidean action is written as 



/[a(r),(/.(T),/3±(T)] 
97r 



4AG 



97r 
2AG 



dr 



C(0,t)) 



W - a(l - V(/3±)) + a302 + ^2^3 ^2 



dra[a2(l + /Lt202)_;L + y(/3±)] 



C(0,v) 



(32) 



Note that the action is given by line integral over the 
contour in complex r-plane, C^.v), start from origin to 

V = X + iY . Following Lyons K|, we integrate along the 
real axis to X first, and then along the imaginary axis to 

V = X + iY where the variables are asymptotic real. 

At the initial position t = 0, we require that the scale 
factor a(0) — and all other variables are regular. Thus 



P2 
M2 



\P2e 



iQr, 



P2e 



iQr, 



\M2\e 



(37) 
(38) 
(39) 



RESULTS 



COMPLEX SOLUTION 



We now solve the complex Einstein's equations with 
the condition that the imaginary part of the dynamical 
variables become zero at large r. 

Figure[T]and Figurel2]show the real and imaginary part 
of the scale factor a, scalar field 0, anisotropy (3±. The 
parameters are /i = 3/4, 0o = 2, P2 = 0.02, and M2 = 
0.15. If P2 = M2 = 0, the parameters are same as HHH's 
model. All imaginary parts become zero at large scale as 
seen by these figures. 

FigureO shows the real part of /3+ vs /?_ in the 
triangular-shaped potential. We can see that the uni- 
verse point bounces between steep walls of the potential. 
However, note that this is an Euclidean solution and not 
directly corresponds to the classical history. In fact, the 
corresponding Lorentzian history has completely differ- 
ent trajectory, see Figure[6l 

The real and imaginary parts of the Euclidean action 
are shown in FigureH) While the real part of the ac- 
tion is oscillating in y < 1, the real part stabilizes and 
imaginary part is rapidly varying in j/ > 2 (in large scale). 
This solution satisfies the classicality condition , thus 
there exist a corresponding classical Lorentzian history. 
The real part of Euclidean action provides the relative 
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probability of history. 

In FigurelH the real part of action stabilizes at Ifi ^ 0.5. 
In the isotropic model, Ir ~ —0.1. We can see that the 
probability of the anisotropic universe is less than that 
of isotropic one. 

B. CLASSICAL LORENTZIAN SOLUTIONS 

Next, we calculate the Lorentzian histories by using 
the complex solution as initial condition. 

Figure [5] shows the two Lorentzian histories. One (we 
call model A) has large anisotropy, the another (model 
B) has small anisotropy. If there is no anisotropy, the his- 
tory shows initial bouncing behavior as shown by HHH. 
Similarly, model B is bouncing as in isotropic case. On 
the other hand, model A has an initial singularity. Usu- 
ally, the existence of singularity is regarded as a sign 
of the breakdown of the theory. However HHH claimed 
that the no-boundary wave function predicts the proba- 
bility of classical histories rather than initial data, thus 
the initial singularity is not a signal of break down of 
the theory. Since the second time derivative of the scale 
factor cPa/dfi is positive in both models, these histories 
have inflationary period. 

The motion of anisotropy in triangular-shaped poten- 
tial is shown in FigurelHl Model A has started with high 
anisotropy, and then rolled down to narrow valley, and 
arrived at low anisotropic state. Model B has started 
at the bottom area of potential, and wandering around 
the initial position. In both models, the anisotropy will 
approach to constant value at large scale, and nearly van- 
ishes. These features are consistent with preceding stud- 
ies (e.g. Hawking and Luttrell [ll|, Wright and Moss (T^ 
etc.). 

The real part of the Euclidean action is shown in 
Figurel?] The figure clearly shows that the probabil- 
ity of high anisotropic histories are smaller than low 
anisotropic histories. Model B's probability is slightly 
smaller than isotropic one. 

FigurelS] shows the area where the classical histories 
exist in P2 — M2 plane. We do not change /i, 4>q. White 
points represent the histories with initial bounce, and 
black points represent histories with initial singularity. 
Clearly we can see that there is the critical line where 
the initial state changed drastically, the high values of 
P2 , M2 area have initial singularity. Note that P2 , M2 
are the absolute values of second differential of /3+, /3_ in 
initial condition, not the initial value. We could not find 
the solution outside of the plotted area. There seems to 
be an upper limit in each P2, M2. 

The relative probability distribution is shown in Fig- 
ure [SI The probability has maximum at isotropic history. 



and is rapidly decreasing with the increase of anisotropy. 
Figure [TU] shows the distribution of /3+, /?- at the bounce, 
and Figured!] shows their probability distribution. We 
can see that states with low anisotropics have higher 
probabilities. 

We found that the distribution of the anisotropy near 
the bounce is highly inhomogeneous. This is very differ- 
ent from classical result. While the classical result pre- 
dicts chaotic behavior near the initial singularity [l^ [3] , 
we could not see such a behavior. 

The anisotropy at large scale are shown in Figure [T^] 
and FiguredHl Clearly, the probability is peaked around 
the isotropic history. Therefore the no-boundary wave 
function of the universe predicts nearly isotropic uni- 
verse. 



VI. CONCLUSION & DISCUSSION 

We investigated the emergence of the classical Mix- 
master universe with a single scalar field and with the 
cosmological constant from the no-boundary quantum 
state according to the method of HHH. We confirmed 
that there are complex solutions that satisfy the classi- 
cality condition, thus there are classical Lorentzian his- 
tories with any amount of anisotropics. This fact im- 
plies that the method of HHH can be applied for a wide 
range of models. It is also found that the probability 
of anisotropic universe is lower than isotropic one. In 
addition the emerged classical universe experiences the 
inflationary expansion and the remaining anisotropy de- 
creases rapidly. Thus the quantum cosmology with no- 
boundary wave function may be able to explain quite 
naturally the emergence of our universe. 

We also found that the classical histories have initial 
singularity or initial bounce depending on the initial ac- 
celeration of anisotropics. High acceleration may induce 
the initial singularity. The no-boundary proposal seems 
to predict that the distribution of anisotropy near the 
bounce or singularity is not homogeneous and very dif- 
ferent from the prediction by classical theory and by anal- 
ysis of Wheeler-DeWitt equation fi\ . On the other hand, 
the loop quantum cosmology predicts non-chaotic behav- 
ior near singularity [lB|- Thus it would be very interest- 
ing to see if the non-chaotic behavior is general in more 
general spacetimes. 
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FIG. 1: The real and imaginary part of the complex-Euclidean solution on X-axis. 
From top to bottom, scale factor a, scalar field (fi, anisotropy /3+ and /3_. 
Left row are real part, right row are imaginary part of the variables. 

H = 3/4,</-o = 2,P2 = 0.02, M2 = 0.15 and 0.871,61 = 0.0942, Op = 3.01,9™ = -0.133 




FIG. 2: The real and imaginary part of the complex-Euclidean solution on Y-axis 
All imaginary parts become zero at large scale. 
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FIG. 3: The real part of /3+ vs 

The universe point comes from high anisotropy to low anisotropy with some bounce. 
Note that this is in Euclidean solution, not classical Lorentzian. 
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FIG. 4: The real and imaginary part of the Euclidean action /. 
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FIG. 5: Lorentzian history of two models. From top to bottom, a, (j), j3+, j3-. 
Left row is model A (high anisotropy), right row is model B (little anisotropy) 
We use the complex Euclidean solution at F = 4 as initial value ai t — 4. 
Model A: fi^ 3/4, 0o = 2, P2 = 0.02, M2 = 0.15 
Model B: /i = 3/4, 0o = 2, P2 = 0.002, A/2 = 0.015 




FIG. 7: The real part of the Euchdean action. Left is model A, right is model B. 
The probability of model A is smaller than model B. 
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FIG. 8: The area where the classical histories exist in P2 — M2 plane. 

White points represents the histories which has initial bounce, black ones has initial singularity. 
There is critical line where the initial state changed and upper limits. 
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FIG. 9: The relative probability distribution in P2 — M2 plane. 

Left represents same area with FigO right represents the area nearby origin (Isotropic). 
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FIG. 10: /3+ vs /3_ at the bounce. 




FIG. 11: The relative probabihty distribution in (3+ — P- plane at the bounce. 
There is high probability surface but no clear structure at low probability. 
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FIG. 12: 13+ vs /3_ at large scale in classical history. 

White region represents bouncing history, black region has initial singularity. 

Since the anisotropy is constant at large scale (see Fig. [5} , this figure represents the observable one. 
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FIG. 13: The relative probability distribution in /3+ — /3_ at large scale in classical history. 
The probability distribution is peaked around the isotropic universe. 



